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Tj- , Abstract 

m 

• ' We apply the ADHMN construction to obtain the SU{n + 1) (for generic values 

in ; 

^^ I of n) spherically symmetric BPS monopoles with minimal symmetry breaking. In 

particular, the problem simplifies by solving the Weyl equation, leading to a set of 

coupled equations, whose solutions are expressed in terms of the Whittaker functions. 

^ ' Next, this construction is generalized for non- commutative SU{n + 1) BPS monopoles, 

C^ '. 

where the corresponding solutions are given in terms of the Heun B functions. 
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1 Introduction 

Bogomolny-Prasad-Soinmerfield (BPS) monopoles are topological solitons in a Yang-Mills- 
Higgs gauge theory in three space dimensions. The equation for static monopoles is in- 
tegrable, and thus a variety of techniques are available for constructing and studying the 
monopole solutions. Direct construction of solutions of the Bogomolny equation with monopole 
number greater than one is a difficult task with the exception of spherical symmetry, (see, for 
example, Ref [1]-[1] and References therein). To bypass this difficulty a number of intriguing 
ideas have been put forward rendering this complicated problem more tractable. 

A powerful approach introduced by Nahm [5j is the so-called Atiyah-Drinfeld-Hitchin- 
Manin-Nahm (ADHMN) construction. The Atiyah-Drinfeld-Hitchin-Manin (ADHM) con- 
struction [6], allows the construction of instantons in terms of linear algebras in a vector 
space, whose dimension is related to the instanton number. Since monopoles correspond to 
infinite action instantons, then an adaptation of the ADHM construction involving an infi- 
nite dimensional vector space might be also possible. Nahm was able to formulate such an 
adaption in the ADHMN construction. To perform this construction a nonlinear ordinary 
differential equation (i.e. the Nahm equation), must be solved and its solutions (i.e. the 
Nahm data), used to define the Weyl equation. 



The Nahin equations provide a system of non-linear ordinary differential equations 

^ = le,MTM,Us)] (1) 

for three n x n anti-hermitian matrices Tj(s) (the so-called Nahni data) of functions of the 
variable s, where n is the magnetic charge of the BPS monopole configuration. The tensor 
Eijk is the totally antisymmetric tensor. These equations may be obtained from the self- 
dual Yang-Mills equations in four dimensions by imposing translation invariance in three 
dimensions. Thus, equations ([T]) are completely integrable, and can be solved in terms of 
abelian functions. 

In the ADHMN approach, the construction of SU{n + 1) monopole solutions to the 
Bogomolny equation with topological charge n is translated to the following problem which 
is known as the inverse Nahm transform |5]. 

Finding the Nahm data effectively solves the nonlinear part of the monopole construction, 
however in order to calculate the fields themselves the linear part of the ADHMN construction 
must be implemented. Given the Nahm data for a n-monopole the one-dimensional Weyl 
equation 

( hn-j In ® Xjaj + iTj ® aj j v(x, s) = (2) 

for the complex 2r2- vector v(x, s), must be solved. I„ denotes the n x n identity matrix and 
X = {xi,X2,X3) is the position in space at which the monopole fields are to be calculated. 
Let us choose an orthonormal basis for these solutions, satisfying 

v"^v ds = I. (3) 

Given v(x, s), the normalized vector computed from ([2]) and ([3]), the Higgs field $ and gauge 
potential Ai are given by 



$ = -i / s\^\ds, (4) 

Ai = f^^diYds, (5) 



where the integrations are to be performed over the range spanned by the minimum and 
maximum eigenvalues of the asymptotic form of $. Then the corresponding gauge and Higgs 



fields satisfy tlie Bogomolny equation and tlie boundary conditions for a charge n monopole, 
and they are smooth functions of x. Although many general results have been obtained few 
explicit solutions are known. Specifically, many solutions of the Bogomolny equation have 
been discovered using the inverse Nahm transform, however this does not mean that analytic 
expressions for the gauge and Higgs fields are known. 

One of the main aims of the present investigation is the construction of charge n, SU{n+l) 
monopoles with symmetry breaking from U{n). This means that the asymptotic value of 
the Higgs field has n equal eigenvalues; which is also referred to as minimal symmetry 
breaking. Moreover, we construct the non-commutative version of our solutions. As was 
recently shown, quantum field theories in non-commutative space-time naturally arise as a 
decoupling limit of the world volume dynamics of D-branes in a constant Neveu-Schwarz- 
Neveu-Schwarz (NS-NS) two-form background [7J. In particular, we explicitly construct the 
solutions of the non- commutative EPS equations to the linear order of the non-commutativity 
scale using a deformed (non-commutative) version of the Nahm equation. 

2 Weyl Equation and SU{n + 1) BPS Monopoles 

In this section, the solutions of the Weyl equation for spherically symmetric SU{n + 1) 
BPS monopoles (for generic values of n) in the minim symmetry breaking case are obtained 
explicitly. By spherical symmetry we simply mean here, that one has to deal only with the 
r dependence, i.e. no azimuthal dependence is considered. To our knowledge the only work 
within this spirit is presented in [HI |9], but again only for particular values of n (i.e. n = 2 
and n = 3, respectively). 

In the minimal symmetry breaking case the Nahm data Tj are defined on a single interval 
[— n, 1] with the only pole occurring at the left-hand end of the interval. This allows a 
construction of Nahm data for charge n monopoles in a minimally broken SU{n + 1) theory 
in terms of rescaled Nahm data for SU{2) monopoles, where the rescaling moves the second 
pole in Nahm data outside the interval. For convenience we shift s so that the Nahm data 
are defined in the interval, [0,n + 1]. Finally, the boundary conditions require that this 
representation in the unique irreducible n- dimensional representation of su{2) (see e.g. [T^ 



for a more detailed discussion). Thus, the Nahm data can be cast as 

Ti = -'-f,n, 2 = 1,2,3 (6) 

where r^'s form the n-dimensional representation of SU{2) and satisfy 

[n, Tj] = 2ieijkTk. (7) 

The n-dimensional representation of SU{2) is of the form 

n— 1 n— 1 n 

A;=l fc=l A:=l 

where e,-" are n x n matrices defined by: ( e,-" ) = 6ik Sji and 
■^ \ ■' J ki 

ak = n + l-2k, Ck = \/k{n-k). (9) 

The Nahm data for SU{n) spherically symmetric monopole of charge n — 1 are given by 
([6]) where fi = f = —- (see also e.g. [HI [10] )■ Since the monopole with this Nahm data 
is spherically symmetric, we may choose {xi,X2,xs) = (0,0, r), and assume that the vector 
v(x, s) is of the form 

n 

v(x, s) = J2 hi{r, s) e!") ® {b,{r, s) ef ^ + h{r, s) ef ) (10) 

1=1 

where e^" is the ra-dimensional column vector with 1 at the position k E Z"*" and elsewhere, 
i.e. the standard basis of M". 

Then the Weyl equation ([2]) becomes 

J, n—l p n—1 

+ ( E c. (4lVi + 4i.) ® (eg + 4?) -{y:c>^ (4i. - 4:V0 » (4? - ^ 
/ 

k=l 1=1 

To proceed with our computation we exploit the following properties: 



A + ZVr. Tp^") +p(") ^6?) ^^(2) , (2)\ _ /V-^ / (n) _ (n) A ^ / J2) _ (2) 

k=l k=l 

J. n n 

+ ^ E -^' 4:^ ® (eS? - eg) - r I ® (ej? - eg) ] ^ /.. eS") (6,eg + h, e^?) = 0. (11) 






With the use of the latter identities and after setting 

ui{r,s) = hi{r,s)bi{r,s), Wi{r,s) = hi{r,s)b2{r,s), (13) 

equation ( ITT1) is equivalent to the following first-order system of differential equations 

ui- ( — ai + r j ui = 0, 

iik+i Ck Wk- { — ctfc+i + r I Uk+i = 0, 

s \2s J 

Wk Ck Uk+i + ( — ctfe + r j Wfc = 0, 

Wn + {:z-an + r)wn = 0, (14) 

where k = 1,2, . . . ,n ~ 1. Here, Ui and Wi for i = 1, . . . ,n are the total derivatives of the 
functions Mj(r, s) and Wi{r, s) with respect to the argument s. Solving these equations is the 
first step in reconstructing a solution of the Bogomolny equations from Nahm data. Then 
the problem of recovering the Higgs and gauge fields is linear. 

Let us now solve these equations. Assume first that r ^ 0. The first and last equations 
may be immediately integrated and their solutions are equal to 

ui = Ki{r)y/s^e'''', Wn = i^2{r) —j^, (15) 



where Kj(r) for i = 1, 2 are the constants of integration. The coupled equations for Uk+i and 
Wk are equivalent (by expressing Uk+i in terms of Wk) to the single second-order equation 

(Tl — 1 \ 
r'^s^ + {n-l-2k)rs^ — jw^ = 0, (16) 

which may be solved by substituting Wk = Wk/s and z = 2rs. The latter equation is then 
reduced to the familiar Whittaker equation: 



1 ^ 

T + - 


2k-n+l) 


1 n2 


2 


■+' / 


4 


Z 





Wk + [-- + ^^ + ^^ ]Wk = Q, (17) 

where W'^ = ^. 

Therefore, the solutions of f lT7|) are given in a closed form, in terms of the Whittaker 
functions as 

Til Tl \ ( Ti \ Ti 

Wk = c^{r)M\ -- + - + /c,-;2rsj + C2(r) l^ (-- + - + /c, -; 2rs 



where Ci{r) for i = 1,2 are constants (see Appendix, for a brief review on Whittaker func- 
tions). In Table 1 and 2 specific examples of the Whittaker M (— | + | + /c, ^; 2rs) and 
W (— I + I + /c, |; 2rs) functions are presented. 



M(-^l+i,t;2..) 


n = 2 


n = 3 


n = 4 


k=l 

k = 2 
k = 3 


2-^/2(smh(rs)-rse-'-=) 


12(rs cosh(rs)— sinh(rs)) 
rs 

6 sinh(rs)-rs(l+rs)e"''= 


12v^ 3smh(rs)-rse-''=+rs(rs-2) e''^ 
(rs)3/2 

12v^ rse''^+rs(rs+2)e"''''-3 sinh{rs) 
(rs)3/2 

4v^3sinh(rs)-rs(3+3rs+2r2s2)e-''=] 
(rs)3/2 



M(2^^^.t;2rs) 


n = 5 


n = 6 


k = 2 
k = 3 
k = A 
k = 5 


20[3rse-''=+rs(9-6rs+2r2s2)e'-s-12sinh(rs) 
(rs)2 

60 6sinh(rs)-rs(rs+3)e"''=+rs(rs-3)e''^ 


60\/2[l5smh(rs)-3rse-''«+rs(r3s3-4r2s2_|_9rs_12)e''= 
(rs)-V2 

60v^ 3rs{rs+4) e-^''+rs(2r'^s'^~9rs+lS) e''=-30sinh(rs)] 
(rs)5/2 

60^2 30sinh{rs)+3rs(rs-4) e'"»-rs(18+9rs+2r2s2) e-''^] 
(rs)5/2 

60v^[3rse''=+rs(12+9rs+4r2s24.j.3s3)e-'-«-15sinh(rs)] 
(rs)5/2 

12^2 15sinh(rs)-rs(l5(rs+l)+5r2s2(rs+2)+2r*s'')e-''=] 
(rs)5/2 


(rs)2 

20 3rse'~''+rs{2r^s'^+6rs+9) e"''=-12sinh(rs)] 
(rs)2 

20[6smh{rs)-rs(3+3rs+2r2s2+r3s3)e-''=] 


(rs)2 



Table 1: Exphcit expressions for the Whittaker function M ( ^^ ^""^\ |; 2rs) for 

n = 2, ...,6 



W'(^^^^.t;2r.) 


n = 2 


n = 3 


n = A 


k=l 

k = 2 
k = 3 


(l+2rs)e"'"^ 
V2rs 


rs 

[l+2rs(l+rs)]e-''= 
rs 


(3+2rs) e-''" 
(2rs)3/2 

[3+2rs(rs+2)] e---' 
V^{rs)3/2 

[3(l+2rs)+2r2s2(3+2rs)] e---' 
v^(rs)3/2 



w[^>'-;^+\:^;2rs) 


n = 5 


n = 6 


k= 1 

k = 2 
k = 3 
k = 4: 
k = 5 


{rsr 


(5+2rs)e-''^ 
(2rs)5/2 

(5+4rs+r2s2)e-'"'' 
V2(rs)5/2 

(l5+18rs+9r2s2+2r-3<)3)e-'-'' 
v^(rs)5/2 

V2(l5+24rs+18r2s2+8r3s3+2r4s4)e-'-= 


(r.)^ 


(rs)5/2 
V2(l5+30rs(l+rs)+20r3s3+10r4s*+4rSs5)e-''^ 


(rs)5/2 



Table 2: Explicit expressions for the Whittaker function W ( ^^ ^""^\ |; 2rs) for 

n = 2 6 



The next step is to derive an orthonormal basis of solutions, and then calculate the Higgs 
field (jlj). An orthogonal basis of solutions is given by 

/ \ / \ 





(vA 


















Vi = 





, Vfc+i = 


Introduce 


the in 


ner product 



Uk+1 

\ / 



< V«;,Va > 



Vn+1 










A; G {!,..., n-1}. (19) 



71+1 



vt v\ ds 



''k^X 



'0 



(20) 



then the required solutions of ([2]) are those which are normalizable with respect to ([3]). It is 
clear that the space of normalizable solutions of ([2]) has dimension n + 1. 

If \i, . . . ,\n+i consist an orthonormal basis for this space then the Higgs field is an 
SU{n + 1) diagonal matrix whose elements are given by 

fW+l 



rn+i 

$K = -W (s - n) vj, v^ (is, K G {!,... ,n + 1}. 
Jo 



(21) 



Recall that, the orthogonal vectors vi and ^n+i are given in terms of the functions ui and 
Wn in (TT5|) . The vectors v^+i for k G {l,...,n— 1} are more complicated and are given in 
terms of the Whittaker functions M (^^^^, f ; 2rs) and W (^^^^, f ; 2rs). In particular, 



the Wk and Mjt+i functions are given explicitly by 



Wk = - Wk, 
s 



Uk+i — yr 

^k 



Wk+( ^^ + r]Wk 



(22) 



where Wk is defined in flTSj) . 

Performing the required integrals fl20|) . and using flTSj) . f l22|) . gives 



^fl 



A4+1 



A/;+i 



n+l 



c! 



^e'"'' c?s, 



1 diy.2 / n-l-2k 

" + r + - 



2 ds 



2s 



W'k 



n+l 
s=0 



he {l,...,n-l}. 



n+l „-rs 



(is. 



(23) 



Note that the normalization factors A/fc, solely consist of boundary terms and thus can be 
easily evaluated. However, the Whittaker W { ^^''~^~^^ , |; 2rs) function at s = tends to 
infinity (see Appendix). Hence setting C2(r) = in (TT5I) to avoid the divergencies, the 
solutions Wk are given in terms of the Whittaker M ( ^^~^""'"^ , ^; 2rs) function only. 



Then the corresponding diagonal elements of the Higgs field at (0, 0, r) are given by 



n+l 



$1 = -—r / {s -n) Vs'^i e""" ds, 
•^1 Jo 



$ 



k+l 



A/fc+1 C*A; 



-Mc+1 C*;. 



sdW, 



k +{rs+'4-k-l]W^ 



n 
2 ds"+l"+2 



n+l 
s=0 



r-n+l 



2s 



r pyf ds + i ra, k e {I, . . . ,n - 1}, 



$ 



n+l 



s — n) 



ds. 



A/'n+l Vo '" '"' Vs^ "^' 

It is a simple task to verify that this solution has indeed the correct asymptotic behavior, 
and also to recover the results obtained for the special cases n = 2 and n = 3 in [SJ [H], 
respectively. 

To conclude, with the method applied in this section starting from the n-dimensional 
representation of the SU{2) algebra (n-monopoles) one ends up with n + l orthonormal 



vectors, which provide the Higgs field corresponding to the minimal symmetry breaking of 
SU{n + 1). Note that we have focused here on spherically symmetric monopoles located at 
(0, 0, r). Similarly, the azimuthal dependence can be implemented in a quite straightforward 
manner via suitable similarity transformations, or by solving the full Weyl equation with the 
simple Nahm data. Then the gauge potentials can also be recovered. This is a mathemati- 
cally and physically interesting problem, and will be addressed in full detail in a forthcoming 
publication. 



3 Non-Commutative SU{n + 1) BPS Monopoles 

The idea of non-commutative space-time offers a smooth way to introduce nonlocality into 
field theories without loosing control. Motivated by string theory the investigation of non- 
Abelian gauge theories defined on non-commutative space-time is of great interest in the last 
few years. In particular, the dynamics of non-Abelian gauge fields involves non-pertrubative 
field configurations (e.g. instantons, monopoles, etc) in an essential way. Thus, in order to 
quantize a gauge theory it is mandatory to study its classical solutions and characterize their 
moduli spaces. 

In [11] Bak derived a deformed Nahm equation for the BPS equation in the non-commutative 
A^ = 4 super-symmetric U{2) Yang-Mills theory. This way, he was able to explicitly con- 
struct a monopole solution of the non-commutative BPS equation to the linear order of the 
non-commutative scale. 

In this section, following Bak, the non-commutative SU{n + 1) BPS spherical symmetric 
monopoles with minimal symmetry breaking, are obtained using a generalization of the 
ADHMN method. The non- commutative construction goes along the same lines as the 
commutative one, with the only difference that all the ordinary products are now replaced 
by the associative (but not commutative) Moyal star product (i^r-product). The latter is 
characterized by a constant positive real parameter 6, which appears in the star product of 
two functions in the following way 



a(x) 7kr6(x) 



gieM-s„a,', 



'" a{x)b{x] 



(25) 



For ^ = it reduces to the ordinary product of functions. In general, 9^^ is a constant, real- 
value antisymmetric {9^^ = —9^^) matrix of order d (where d is the dimension of space-time) . 
For simplicity, we choose its only non-zero terms to be 9i2 = —6*21 = 9. 

In the non-commutative case, the derivation goes through once all the product operations 
are replaced by ^-product. Namely, the Weyl equation ([2]) is modified as 



I — - — h I„ ® Xjaj — iTj (g) aj j ^ v(x, s) = 0. 



(26) 



Ti satisfy (for a more detailed discussion on these issues see [TT]) the deformed Nahm equa- 
tion: 

-^ = 2^ufcK, Tk]-96,3l. (27) 

Setting Ti = Ti — 9s6i3, is easy to see that Tj satisfy the Nahm equations ([1]); they have 
simple poles; and their residues are irreducible representations of SU{2). 
The deformed Weyl equation takes then the following form 

( -— + In ® Xjaj - ifj ®aj + 9d] v(x, s) = 0, (28) 



where 



= ^In(E)i(Tid2 -cT2di)+ is In ^ as. (29) 



Next we construct the actual solutions of the non- commutative BPS equations from the 
deformed Nahm data Ti fl27|) . As before, we focus on the spherically symmetric case by 
setting Xi = (0, 0, r), and Tj are the ordinary Nahm data of the form ([6]) where fi = f = — -• 
In this case, the deformed Weyl equation ( l28l) becomes 

{-—-^-Ti®ai + fl®aAw{f,s) = Q, (30) 

where f = r + i9s. 

Similarly to the commutative case, the deformed Weyl equation f l30|) leads to a first-order 
system of differential equations given by f lT^ for r —> r. Then the first and last equations 
are decoupled and their solutions are given by 

if>^2\ , ._. e V ^ 2 



Ml = A;i(f)Vs"i e^"'+-' >, Wn = k2{f) ^^. (31) 

10 



Finally, the coupled equations for Uk+i and Wk are equivalent to the deformed single second- 
oder differential equation 



s idk + 2swk —Irs— iOs + {n — 2k — l)fs + 



Wk = 0, 



(32) 



for k & {1, . . . ,n — 1}. Note that, the aforementioned equations become the corresponding 
equations of the previous section when 6 = 0. 

The latter equation is solved in terms of the so-called Heun B functions denoted hence- 
forth as Hb (see. Appendix, for more information regarding Hb and the corresponding 
differential equation). The solution is given by 

--,.,,, , 2(-l)ir , (-l)i(4A; + 2-2n)r , ,3 /- 



Wk = Ci(r)v^s"'~^ e^^^ 






^'^"^'^^Hb I -n. 



2(-l)3r 
V~9 ' 



n-2A;-2, 



l)i(4^+2-2n)r 3^; 



v^ 



(33) 



where (— 1)* = H= and (—1)* = —^- It can be shown (see Appendix) that for the specific 
form of the solutions given by f l5^ the corresponding Heun B functions satisfy the conditions 



Hb i-n) 



;-i)tv^5 



Hnin). 



(34) 



Therefore, the general solution (13311 becomes 



.-1 .^rs^^s^) u. i . '±}}IL, . _ 2. _ 2, (-1)^(4^ + 2 -2.)r^ ^_^^ 3 ^^ 



Wk = c{r,e)Vs^~^ey''^'^'')HB I n. 



Ve 



y/e 



(35) 



where c{r, 6) is an arbitrary function of r and 6. Table 3 presents the series expansion of 
the Heun B function Hb [n 
order: 0{s^). 



2ti^,„-2-2*, 



(-l)./.(4|+2-2„).- . (_l)S/4^,A „p t„ third 



HB{n) oc O (s3) 


n = 2 


n = 3 


n = 4 


k = l 
k = 2 
k = 3 


1 - |rs + (r2 - |/0) s2 


1 - rs + f (r2 - 70) s2 
1 _ l^s + I (3^2 _ 210) s^ 


l-lrs+[y-lie)s^ 



11 



Table 3: The Hb (n, 



order Ofs^) for n = 2 4. 



l)^/^V^s) function up to 



Next in order to obtain the deformed Higgs field we follow the procedure of the pre- 
vious section. First we choose the deformed orthogonal basis similarly to ( IT9|) . Then the 
normalization factors are equal to 



M 



A4+1 



s"i e 



/^+f ^^) ds. 



CI 



^ + {r + ies + z I Wf 



Afn 



+1 



2 ds 



2s 



n+l 
s=0 



ke {l,...,n-l}, 



ds. 



(36) 



'0 ys"" 

where VF^ = swk and w^ is given by fl35|) . The corresponding diagonal elements of the Higgs 
field are then expressed as 

cn+l 



$1 = - 



M 



fs-nlA/s^e*-""^^'') rfs, 



$fc 



fc+i 



■Mc+iC'fc 



st/Vr,2 



(77/ 
rs + z6's^ H A; - 1 



A/fc+iC*^ Jo 

r-n+l 



2 rfs 

2^ 



W^.^ 



n+l 
s=0 



r -i9s]W^ ds + in, k e {I,. . .,n- 1}, 



$n+l 



(s — n) : — (is. 



(37) 



A^n+l Jo ' ' VS"" 

As expected, the expressions above reduce to equations fl23l) and flM|) when the deformation 
parameter vanishes, i.e. in the limit ^ — t- 0. 



4 Discussion 

We have been able to derive, using a quite elegant and unifying methodology, explicit 
expressions for spherically symmetric commutative and non-commutative SU{n + 1) BPS 
monopoles in the minimal symmetry breaking case. The use of the generic n-dimensional 
representation of SU{2) allowed the explicit construction of solutions of the Weyl equation 



12 



giving rise to sets of simple differential equations. The solutions of the aforementioned dif- 
ferential equations were expressed in terms of confluent (Whittaker) or biconfiuent (Heun B) 
hypergeomentric functions for the commutative and the non- commutative case, respectively. 
Although we have restricted our investigation in the spherically symmetric case it is im- 
portant to note that azimuthal dependence may be implemented in the generic situation 
by means of suitable similarity transformations, however this issue will be discussed in full 
detail in a forthcoming work. It is also worth noting that the case n — )■ oo merits spe- 
cial investigation, given that the behavior of the discovered solutions in this situation is 
rather particular, modifying the form of the final expressions for the Higgs field. Physical 
interpretation of the behavior of the Higgs field is also desirable. Such results may be also 
of relevance when one considers, within this context, infinite dimensional representations of 
SL{2, R) (the non-compact case, e.g. non-trivial spin s = representation of SL{2, R)). The 
solution of the Weyl equation in this case involves a particular polynomial basis, given that 
the associated representation is expressed in terms of simple differential operators. These 
are quite intriguing issues, and will be left for future investigations. 
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A Whittaker and Heun B functions 

The Whittaker functions M{k, m; z) and W{k, m; z) solve the Whittaker differential equation 
(see e.g. [El [13]): 

where 

W = ciM{k,m]z)^C2W {k,m]z). (39) 

They can be defined in terms of confluent hypergeometric functions as follows 

M{k,m,z) = e~^z"'+^iFi(-+m-k;l + 2m;zj, (40) 

13 



W{k,m,z) = e "^ z'' 2F0 i- + m- k,- -m- k;-z '^j , (41) 

where iFi (a; b; z) = e^ iFi {h — a; h; —z) and nFm are generalized hypergeometric functions 
(see e.g. [HI [H [S]) defined as: 



00 



^ {bi)i{b2)i...{bg)i n r(a) 

In the special case where m + ^ — k is a positive integer (which is our case), the Whittaker 
W {k, m, z) function is defined by the integral function 

W (k, m- z) = ,, / r^- 2+-(i + 1)'=- 2+-e-* dt (43) 

r(i-fc + m)Jo ^ 

for all values of z except negative real values. Then, it is obvious that at 2; = (corresponding 

to s = in our case) it diverges. 

Let us now introduce the Heun Biconfluent differential equation 

d^y ( \^a\dy ( 1 (a + 1)6 + d^ 



,.H-^-2--6+— j^+^.-a-2--^.^Jy = 0. (44) 

The Heun B function Hsia, b, c, d; z) is a Frobenius solution to Heun's Biconfluent equation 
(jUj) computed as a power series expansion around the origin, which is a regular singular 
point. In this case, the growth of solutions is bounded, in any small sector, by an algebraic 
function. Because the next singularity is located at infinity, this series converges in the whole 
complex plane. 

Note that the Heun B function can be expressed in terms of the Whittaker M function 
for specific values of the parameters, i.e. 

HB{a, 0, c, 0; 2;) = -^ e^ M (|, ^; z^^ . (45) 

It is thus clear that the results of section 3 reduce to those of section 2, in the limit — ;■ 0. 

A special case occurs when in Hsia, b, c, d, z), the third parameter satisfies the condition 

c = 2(k + 1) + a where k is a positive integer (which is our case). In this case the k^ + 1 

coefficient in the series expansion is a polynomial of degree k in d. When d is a root of this 

th 

polynomial, the k^ +1 and subsequent coefficients cancel and the series truncates resulting 
in a polynomial form of degree k for Heun B. 
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